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The  application  of  the  eigensystem  realization  algorithm  to  flexible  structures  excited  with  nonimpulse 
type  inputs  is  described.  Sinusoidal  pulses  are  a  special  subset  of  such  inputs,  and  a  discussion  is  given 
on  their  use  to  lock  on  to  particular  modes.  The  determination  of  impulse  response  functions  from 
nonimpulse  response  functions  is  performed  via  a  two  step  procedure.  In  the  first  step,  free  decay  data 
is  analyzed  by  the  eigensystem  realization  algorithm  to  determine  the  structural  frequencies  and  damping. 
Forced  response  data  is  analyzed  in  the  second  step,  by  solving  an  overdetermined  set  of  linear  equations 
for  the  modal  coefficients  of  the  impulse  responses.  It  is  shown  that  the  number  of  impulse  response 
coefficients  is  halved  if  displacement,  velocity,  and/or  acceleration  sensors  are  used  in  conjunction  with 
force  actuators.  System  identification  is  completed  after  the  eigensystem  realization  algorithm  is  used  on 
the  impulse  responses  to  produce  a  linear  state  space  model  of  the  structure.  Numerical  examples  are 
given  in  which  structural  models  are  determined  from  simulated  data  corrupted  with  noise. 


Introduction 

HE  eigensystem  realization  algorithm  (ERA)  is  a  system 
identification  method  that  has  been  studied  extensively.1 ' 
It  determines  a  linear  differential  equation  model  of  a  structure 
by  analyzing  the  structural  impulse  response  data.  The  applica¬ 
tion  of  an  impulse  to  a  structure  can  be  a  problem,  since  it 
excites  all  of  the  structural  modes.  A  model  that  is  identified 
from  this  data  would  be  larger  than  necessary,  if  the  structure 
were  to  be  operated  in  a  band-limited  manner. 

In  practice,  one  might  apply  a  short  pulse  to  the  structure  that 
would  approximate  an  impulse  for  the  low-frequency  modes  and 
would  not  excite  the  high-frequency  ones.  This  procedure  would 
be  successful  if  the  magnitude  of  the  pulse  were  enough  to 
produce  a  significant  signal  at  the  sensor  outputs.  If  this  is  not 
the  case,  then  a  more  frequency  band-limited  signal  would  be 
a  logical  choice  of  input.  For  example,  several  cycles  of  a 
sinusoidal  wave  at  or  near  one  of  the  resonant  frequencies  could 
be  applied.  A  lightly  damped  structure  may  not  survive  extended 
excitation  at  one  of  the  resonant  frequencies;  the  sinusoidal 
pulse  should  be  stopped  after  sufficient  excitation  has  been 
achieved.  By  conducting  a  set  of  experiments  involving  sinusoi¬ 
dal  pulses  of  different  frequencies,  the  structural  modes  within  a 
certain  frequency  range  can  be  observed.  Finite  element  analysis 
and  previous  system  identifications  can  be  used  to  narrow  dow n 
the  frequency  bands  of  interest.  ERA  cannot  be  used  to  analyze 
this  nonimpulse  response  data;  an  algorithm  is  needed  to  deter¬ 
mine  the  impulse  response  from  nonimpulse  response  data.  One 
such  algorithm  was  developed  in  recent  work  by  Horta  et  al.,3 
Phan  et  al.,4  Juang  et  al.,5  and  Phan  et  al.6  The  algorithm,  known 


as  Observer/Kalman  Filter  Identification,  identifies  an  observer/ 
Kalman  filter  to  determine  the  impulse  response. 

The  remainder  of  this  paper  describes  an  alternate  method 
that  does  not  involve  an  observer/Kalman  filter.  Instead,  ERA 
is  first  used  to  determine  the  frequencies  and  damping  from 
the  unforced  response.  The  impulse  response  can  then  be  found 
from  the  least  squares  solution  of  ar.  overdetermined  set  of 
linear  equations.  The  impulse  responses  between  each  input 
and  output  pair  are  then  simultaneously  analyzed  by  ERA  to 
produce  a  linear  model  of  the  structure. 

Derivation  of  Procedure 

Consider  the  following  model  of  a  flexible  structure: 

Mq  r  Cq  +  Kq  =  flu 

„  (1) 
yd  =  Ctq  ,  y,  =  C,q  ,  y„  =  Caq 

where  the  displacement  q,  the  input  force  u,  the  displacement 
sensor  output  yd,  the  velocity  output  y„  and  the  acceler-ation 
output  y0  are  dimensioned  as:  qeR\  us.  Rm,  yd e Rp,  y, e Rp, 
y«  e  Rp. 

Without  loss  of  generality,  we  will  consider  each  input-output 
pair  separately  (m  =  1 ,  p  =  1 ).  The  single  output  is  assumed 
to  be  a  displacement,  a  velocity,  or  an  acceleration.  The  mass 
matrix  M  is  assumed  to  be  positive  definite  and  symmetric;  the 
damping  matrix  C  and  the  stiffness  matrix  K  are  assumed  to 
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be  positive  semidefinite  and  symmetric.  It  is  assumed  that  the 
system  is  modalizable, 

v  +  C'i>  +  K'v  =  |3'u 

yj  =  O .  y.  =  cy,  y„  =  civ 

The  matrices  C'  and  K'  are  diagonal.  For  undamped  systems 
and  for  systems  with  Rayleigh  damping,  modalization  can  be 
performed.  This  modalizability  assumption  is  not  essential  for 
the  system  identification  method  described  subsequently,  but 
it  will  be  seen  that  the  number  of  parameters  which  need  to  be 
determined  is  reduced  if  modalization  occurs. 

Consider  one  of  the  modes  and  the  contribution  to  the  output 
from  this  mode, 

V,  +  c'v,  +  Jfc'Vj  =  P'u 

* 

yj,  =  civ,,  yd  =  y  y* 


C  =  (0  c'v).  For  an  acceleration  sensor,  C,  =  {-c',k[  - 
c^c'i).  The  system  described  by  Eq.  (5)  is  the  diagonalized 
version  of  the  system  described  by  Eq.  (6).  Therefore,  the 
transfer  functions  are  equal, 

C,,(sl  -  =  Cisl  -  A.y'B,  (7) 

After  some  algebra,  the  Eq.  (7)  yields  conditions  for  the  dis¬ 
placement  and  the  velocity  sensors. 

Displacement  sensor: 

c,  =  0  (8) 


Velocity  sensor: 


w,c,:  =  -a  c„ 


Acceleration  sensor: 


>V,  =  c,v,, 


><.,  =  CaVi 


=  -  clklvi  +  c'.p/w 


A  i^i  i 

=  2c  + 

*  c,!  c j\  +  ml 


where  the  constants  c,„  c,„  crh,  and  u>„  were  defined  in  Eq.  (4), 
and  the  new  constants  /„  and  /„  are  defined  as  follows: 


Since  it  is  possible  to  perform  the  transformation  v,  =  and 
then  cancel  P'  from  the  given  differential  equation,  we  will 
assume  from  this  point  onward  that  this  has  been  done  by 
setting  P'  equal  to  1 . 

The  impulse  response  from  this  second-order  system  is 
given  by 

Kff)  =  c,e»  +  C(V*'  +  c'b(r) 
s,  =  cj,,  +  jm,r  c,  =  c,,  +  jch 

where  the  asterisk  is  the  complex  conjugate  operator;  c,  and  s, 
are  the  complex  modal  coefficient  and  the  continuous  time 
eigenvalue,  respectively;  and  5(0  is  the  Dirac  delta  function. 
A  first-order  state  space  system  that  has  this  impulse  response 
is  given  by 

Ji,,(0  =  A,xSi(t)  +  B,u(t) 
y,(t )  =  Cltx,,(t)  +  D,u(t) 

*-(:  -)•  -(i)  <5> 

Q  =  (c,  cf),  D,t  =  cl 

The  system  described  by  Eq.  (3)  can  also  be  placed  in  a 
first-order  state  space  form 

x, (t)  =  A,x,{t)  +  B,u(t) 

y, (0  =  C,x,(t)  +  D,«(0 

',i  ■  (-*■  D,'ci 
For  a  displacement  sensor,  C,  =  (cj,  0).  For  a  velocity  sensor, 


/,,  =  1  -  [2<t JAai  +  0)5,)] 

4  =  [2crs,0)s,/(cr5,  +  oj;,)] 


Now  consider  the  discrete  time  version  of  the  system 
described  by  Eq.  (5) 

xr.(k  +  1)  =  F,xr{k)  +  C,M(k) 

y,(k)  =  Hr.xr,(k)  +  (14) 


0  if  ’ 


H„  =  (c,  cf),  <t>,  =  cl 

r,  =  ^  (15) 

=  cr„  +  jmri  (16) 

=  (17) 

where  h  is  the  time  step  size.  If  an  impulse  is  applied  at  k  = 
0  to  this  discrete  time  system,  then  the  discrete  time  impulse 
response  is  equal  to  the  impulse  response  of  the  continuous 
time  system  (5)  at  the  sampling  points: 

hr.(k)  =  Cir!  +  cf(rfy  (18) 

=  cl  (k  =  0)  (19) 

This  is  also  the  form  of  the  unforced  response  (nonzero  initial 
conditions  and  zero  input).  Therefore,  an  analysis  of  the 
unforced  response  by  ERA  will  correctly  determine  the  discrete 
time  eigenvalues.  The  eigenvalues  of  the  state  space  model 
produced  by  ERA  are  the  discrete  time  eigenvalues  r,.  By 
performing  a  complex  logarithm  on  Eq.  (15),  we  can  obtain 
the  continuous  time  eigenvalues  from  the  discrete  time  ones. 
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After  substituting  for  c,  [Eq.  (4)1  and  rt  [Eq.  (17)]  into  the 
equation  (18)  and  simplifying,  we  obtain 

h.,(k)  =  c,ty„(k)  +  c„yn(k)  (20) 

=  <  (k  =  0) 

-y,  i  (Jt)  =  2m,*cosOt0,)  (21) 

ya(k)  =  -2m*  sin(M),) 

Using  conditions  (8),  (9),  and  (10),  we  can  simplify  the 
equation  (20)  for  k  >  0  if  the  modal  frequency  u>„,  is  not  equal 
to  zero,  and  if  /„  is  not  equal  to  zero. 

Displacement  sensor: 

=  cd-fatt)]  (22) 

Velocity  sensor: 

=  c,,  [y,,(k)  -  (o-,/<i>,,)  yn(*)]  (23) 

Acceleration  sensor: 


Real  Matrix  Representation 

A  state  variable  transformation  can  be  used  to  produce  a 
system  representation  of  Eq.  (14),  which  is  composed  of  real 
matrices 


xr.(k)  «  Tx'r.  (k) 

(30) 

H-i !) 

(31) 

f;,  =  T-'F,.T 

(32) 

_  (  <T„ 

<TrJ 

(33) 

G,'  =  T-'Gn 

(34) 

■0 

(35) 

=  h„t 

(36) 

=  (~  2cq  2c, i) 

(37) 

V*)  -  c*  [MA>  7.i (*)  +  7 a(*)]  (24)  Numerical  Examples 

Single  Input-Single  Output 

The  displacement  output  is  equal  to  the  convolution  of  the  The  identification  of  two  structural  modes  is  simulated  for 

impulse  response  (from  all  of  the  modes)  with  the  input  force,  a  displacement  sensor.  From  the  discussion  in  the  previous 

provided  there  are  zero  initial  conditions  section  and  the  displacement  sensor  condition  (8),  the  following 

is  an  example  of  a  state  space  representation  for  a  single  input- 
yd(k)  =  hr(k)  *  u(k)  ■  h  (25)  single  output  structure: 


[  J  v*)|  *  “(*) '  h 

(26) 

/  0.960 

r  f  -0.282 

0.282 

0.960 

0 

0 

0  \ 
0 

*-i  -  1  J 

0 

0.978 

0.208 

it 

\  0 

0 

-0.208 

0.978 

V  cl2[y,2(k)  *  u(k)  ■  h) 

(27) 

\ 

/ 

i*  i 


where  •  h  indicates  multiplication  by  the  time  step  size. 

A  system  identification  experiment  may  run  from  k  =  0  to 
k  =  L.  If  the  number  of  data  points  (L)  is  laiger  than  the  number 
of  modes  (n)  then  Eq.  (27)  is  an  overdetermined  set  of  linear 
equations  which  is  to  be  solved  for  the  coefficients  ca  (i  =  1, 
. . . ,  n).  The  QR  factorization  can  be  used  to  accurately  obtain 
a  least  squares  solution  to  an  overdetermined  set  of  linear 
equations,  as  described  in  Ref.  7. 

The  response  for  a  velocity  sensor  (for  zero  initial  conditions) 
can  be  derived  similarly  and  is  given  by 


o.-  T 

\  0.002/ 

tf,  =  ( 7  0  7  0)  O,  =  0 

As  the  F|  matrix  shows,  the  discrete  time  eigenvalues  are  at 
0.960  ±  j 0.282  and  at  0.978  ±  y0.208.  Any  function  could 
have  been  applied  to  the  input,  but  for  the  purpose  of  this 
example  a  sinusoidal  pulse  was  applied  (Fig.  1).  The  displace- 


*  u(k)  ■  h 


(28) 


with  c,i  (i  =  1, ....  n)  being  the  coefficients  to  be  solved  for. 

Likewise,  the  response  for  an  acceleration  sensor  (for  zero 
initial  conditions)  is  given  by 


-  i 


(k)  +  y„(k) 


After  the  cn  or  ca  coefficients  have  been  determined  for  each 
input-output  pair,  the  impulse  response  for  each  input-output 
pair  can  be  computed  for  k  =  0, ....  L  by  using  Eqs.  (22-24) 
and  summing  over  all  of  the  modes.  The  eigensystem  realization 
algorithm  can  then  be  used  on  all  of  the  impulse  responses 
simultaneously,  to  yield  a  linear  state  space  model  of  the  flexible 
structure  (see  Ref.  1  for  details). 
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ment  output  with  5%  added  noise  is  shown  in  Fig.  2.  Since  the 
input  is  zero  for  k  >  50,  the  eigensystem  realization  algorithm 
was  applied  to  the  unforced  response  in  Fig.  2  for  it  >  50  to 
determine  the  modal  damping  and  frequencies.  The  impulse 
response  coefficients  in  Eq.  (27)  were  solved  for,  and  the 
resulting  modalized  state  space  system  is  given  by 


is  applied  to  the  first  input  (Fig.  3),  with  the  second  input  being 
zero.  The  displacement  outputs  with  10%  added  noise  are  shown 
in  Fig.  4.  In  the  second  experiment,  a  sinusoidal  pulse  is  applied 
to  the  second  input  (Fig.  3),  with  the  first  input  being  zero. 
The  displacement  outputs  with  10%  added  noise  are  shown  in 
Fig.  5. 


0.95  .>7  0  0 

-0.28,  u.959  0  0 

0  0  0.974  0.204 

0  0  -0.204  0.974 


_  _  0.002 
G2  -  o 


H2  =  (6.49  0  8.17  0)  <t>2  =  0 


The  worst  error  occurs  in  the  comparison  b  ,K'  third 

elements  of  the  H,  matrix  and  H,  matrix;  th.  is  \ 


Two  Inputs-Two  Outputs 

As  another  example,  consider  the  case  of  two  structural 
modes  with  identical  frequencies.  The  following  is  a  discrete 
time  system  for  two  displacement  sensors  and  two  force  actua¬ 
tors: 


/  0.951  0.309  0  0 

...  _  -0.309  0.951  0  0 

~  0  0  0.951  0.309 

\  0  0  -0.309  0.951 


0  0 

„  _  2  X  1(T3  10~3 

“  0  0 
,  5  X  10-*  2  X  10  3 


Input  to  One  Actuator 


7  0  3  0 
5  0  7  0 


As  the  F\  matrix  shows,  the  repeated  discrete  time  eigenvalues 
are  at  0.95 1  ±  y0-309.  In  the  first  experiment,  a  sinusoidal  pulse 


10  20  30  40  30  60  70 

Discrete  Tune 

Fig.  3  Input  excitation  at  one  actuator. 


Response  to  Input  at  Fust  Actuator 


10  20  30  40  50  60  70  80  90 

Discrete  Time 

Fig.  4  Noisy  outputs  due  to  first  actuator. 


Single  Input  -  Single  Output 


Response  to  Input  at  Second  Actuator 


-0.4; 
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Since  the  input  is  zero  for  k  >  4 1 ,  the  eigensy stem  realization 
algorithm  was  applied  to  the  unforced  response  in  Fig.  4  and 
5  for  k  >  4 1  to  determine  the  modal  damping  and  frequencies. 
The  impulse  response  coefficients  in  Eq.  (27)  were  solved  for. 
With  these  coefficients,  the  four  impulse  response  functions 
between  the  two  inputs  and  two  outputs  were  computed  using 
Eq.  (22).  The  eigensystem  realization  algorithm  was  then  used 
on  the  four  impulse  response  functions  to  obtain  a  discrete  time 
state  space  model.  After  a  modalizing  transformation  of  that 
state  space  model,  we  obtain 


(  0.950  0.308 
-0.308  0.950 
0  0 

0  0 


0  0 

0  0 

0.950  0.308 
-0.308  0.950 


(4.13  X  lO'" 
2.11  X  10'3 
-4.84  X  lO'" 
5.27  X  10-" 


4.75  X  10" 
1.05  X  lO’3 
-5.40  X  10-" 
2.11  X  lO'3 


/7.00  -5.50  X  10-' 
\5.00  5.26  X  lO"* 


3.00  -5.35  X  10'*\ 
7.00  4.03  X  10"* j 


<t2  = 


In  comparing  the  parameters  of  the  first  model  with  the 
identified  model,  we  find  that  the  error  is  less  than  6%. 

Discussion 

System  identification  using  nonimpulsive  inputs  has  been 
described  in  this  paper.  Band-limited  white  noise  could  be 
used  as  the  nonimpulsive  input,  but  the  sinusoidal  pulse  was 
considered  because  of  its  simplicity.  An  initial  frequency  range 
for  the  sinusoidal  pulse  frequency  can  be  determined  via  finite 
element  analysis  and  previous  system  identifications.  The  fre¬ 
quency  range  is  divided  into  a  grid  with  the  desired  frequency 
resolution,  and  a  set  of  system  identifications  is  performed  for 
each  grid  point.  Any  mode  that  is  detected  by  one  of  these 
system  identifications  can  be  identified  once  more,  by  setting 
the  frequency  of  the  sinusoidal  pulse  to  the  identified  mode 
frequency.  This  causes  the  sinusoidal  pulse  to  "lock  on"  to  the 
detected  mode,  thereby  maximizing  the  signal-to-noise  ratio 
for  that  particular  mode.  The  system  identified  model  increases 
by  at  least  one  mode  whenever  a  lock-on  condition  occurs. 

Passive  vibration  control  is  well  suited  for  use  with  this  form 
of  modal  system  identification.  Analytical  formulas  have  been 
derived  for  passive  vibration  suppression  by  an  inertial  reaction 
device,*  or  by  piezoelectric  actuators.9  After  the  sinusoidal  pulse 
excitation,  the  modes  at  or  near  that  frequency  can  be  identified 
via  the  procedure  in  this  paper.  With  this  identified  model,  the 
analytical  formulas  can  be  used  to  determine  a  passive  vibration 
controller.  This  ease  of  design  and  the  unconditional  stability 
of  passive  controllers  provide  a  means  to  suppress  the  modes 
that  are  excited  by  a  sinusoidal  pulse  at  a  particular  frequency. 
Rapid  suppression  is  desirable,  so  that  the  system  identification 
can  quickly  proceed  to  the  next  frequency  data  point. 

If  modalizability  of  the  second-order  structural  differential 
equations  can  be  assumed,  the  number  of  impulse  response 
coefficients  to  be  determined  is  reduced  by  one-half.  This 
removal  of  parameters  ensures  that  certain  physical  constraints 
are  satisfied.  For  example,  a  displacement  sensor  on  a  flexible 
structure  with  zero  initial  energy  must  show  zero  displacement 
at  time  zero,  even  when  a  force  impulse  is  applied  at  time  zero. 
The  first  modal  coefficient  cf|  is  absent  in  the  modal  component 
of  the  displacement  impulse  response  [Eq.  (22)];  the  presence 
of  the  sine  function  in  the  equation  causes  the  modal  displace¬ 
ment  to  be  zero  at  time  zero.  Since  the  displacement  impulse 


response  is  the  sum  of  the  modal  displacement  impulse 
responses,  it  is  guaranteed  to  be  zero  at  time  zero.  The  situation 
for  velocity  or  acceleration  sensors  is  somewhat  similar,  in  that 
only  one  of  the  two  modal  coefficients  is  to  be  solved  for. 

In  the  case  where  the  damping  is  large  and  modalizability 
does  not  occur,  there  is  no  relationship  between  the  two  modal 
coefficients  for  each  mode.  The  simplifications  presented  in 
the  given  analysis  are  not  performed  in  this  case,  and  both 
modal  coefficients  are  solved  for.  It  is  still  possible,  however, 
to  constrain  the  displacement  impulse  response  to  be  zero  at 
time  zero.  Setting  the  discrete  time  k  to  zero  in  the  expression 
for  the  displacement  impulse  response  will  result  in  a  linear 
equation  that  constrains  the  modal  coefficients.  This  equation 
is  to  be  appended  to  the  overdetermined  set  of  linear  equations 
described  in  this  paper,  and  the  whole  set  of  equations  is  to  be 
solved  for  the  modal  coefficients.  In  the  case  of  a  velocity 
sensor,  the  impulse  response  expression  contains  velocity  modal 
coefficients.  Integrating  the  expression  results  in  a  displacement 
impulse  response,  which  is  constrained  to  be  zero  at  time  zero. 
Thus  a  linear  constraint  is  established  between  the  velocity 
modal  coefficients;  this  is  appended  to  the  overdetermined  set 
of  equations.  A  linear  constraint  is  similarly  obtained  in  the 
acceleration  sensor  case  by  performing  two  integrations  on  the 
impulse  response  expression. 

Structural  models  were  considered  in  this  paper;  the  more 
general  case  of  a  set  of  first-order  differential  equations  can  be 
considered  as  well.  In  this  case,  terms  of  the  form  c.rf  appear 
in  the  discrete  time  impulse  response,  where  r,  is  the  ith  discrete 
time  eigenvalue  and  c,  is  a  coefficient.  However,  terms  of  the 
form  c,kJr*  are  necessary  if  Q  repeated  eigenvalues  occur  (0 
~  j  <  Q)-  The  combination  of  the  eigensystem  realization 
algorithm  and  a  least  squares  determination  of  the  impulse 
response  coefficients  would  also  determine  a  linear  state  space 
model,  for  the  case  of  first-order  differential  equations. 

For  multi-input/multi-output  systems,  the  proposed  proce¬ 
dure  requires  one  input  excitation  at  a  time  and  that  experiments 
be  done  separately  for  each  input-output  pair.  Simultaneous  use 
of  all  inputs  can  be  advantageous  in  that  only  one  experiment 
need  be  performed  for  all  inputs  and  outputs.  However,  one 
experiment  provides  less  data,  which  would  tend  to  increase 
the  sensitivity  of  the  system  identified  model  to  noise.  If  it 
is  critical  to  perform  only  one  experiment,  then  the  system 
identification  procedure  would  need  to  be  slightly  modified. 
The  impulse  response  expressions  would  be  expanded  to  sum 
the  modal  components  from  each  input,  and  more  modal  coeffi¬ 
cients  would  be  simultaneously  solved  for. 

A  more  significant  reason  for  the  simultaneous  use  of  all 
inputs  is  to  increase  the  excitation  of  the  structure.  In  the  case 
of  many  weak  actuators  distributed  over  a  large  structure,  the 
use  of  only  one  actuator  may  not  produce  significant  sensor 
signals.  An  input  transformation  can  be  performed  where  new 
inputs  are  defined  to  be  linear  combinations  of  the  physical 
inputs.  Experiments  would  use  only  one  of  the  new  inputs  at 
a  time,  but  this  would  imply  the  simultaneous  use  of  all  of  the 
physical  inputs.  With  the  significant  sensor  signals  resulting 
from  the  simultaneous  input  usage,  a  model  is  identified  that 
relates  the  new  set  of  inputs  to  the  sensor  outputs.  A  model 
relating  the  physical  inputs  to  the  sensor  outputs  can  be  obtained 
by  inverting  the  input  transformation. 

Concluding  Remarks 

A  method  has  been  presented  that  applies  the  eigensystem 
realization  algorithm  to  flexible  structures  excited  with  nonim¬ 
pulse  type  inputs.  The  simplicity  of  the  method  makes  it  an 
interesting  alternative  to  the  OKID  method.  However,  it  is  more 
complex  than  the  determination  of  the  impulse  response  via  the 
inverse  Fourier  transform  of  the  frequency  response  function. 

Both  steps  of  the  presented  method  possess  least  squares 
characteristics  that  minimize  the  sensitivity  of  the  identified 
model  to  noise.  The  inverse  Fourier  transform  method,  with  its 
simplicity,  does  not  attempt  to  minimize  noise  sensitivity.  A 
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quantitative  comparison  of  these  three  methods  with  respect  to 
simplicity  and  noise  sensitivity,  therefore,  seems  appropriate 
for  future  research. 
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